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ABSTRACT 


This work is a study of processing techniques for 
presenting continuous tone or computer generated multi gray 
level pictures on bilevel displays. These displays are 
bilevel in nature and have individual display cells, all 
of the same size, arranged in a rectangular array. Six 
processing techniques have been discussed. All the 
techniques achieve the subjective effect of continuous 
tone by properly controlling only the spatial density of 
bilevel display states, 

A nonlinear model of the human visual system has 
also been considered to show analytically the eye's 
perception of spatially arranged binary patterns as 
different gray levels. Two sets of such binary patterns 
have then been used to display a computer generated 
picture. 

In the remaining processing techniques the input 

image is divided into picture elements (pixels) and 
intensity of each pixel is then compared with a threshold 
value. The corresponding- output pixel is turned to 
bright state if the intensity is greater than the threshold. 
The processing techniques differ In the manner the threshold 
value is calculat^. Images processed by all the techniques 
are exhibitedj^^qv^uated compared. 


CHAPTER I 


INTRCBUCTION 

This work is a study of techniques/ algorithms that can 
be used to transform a multi-level continuous tone/computer 
generated still, monochrome image for reproduction on 
binary displays/printing systems. In other words, the 
algorithms transform a multi-gray level picture into a 
bilevel picture which retains the same gray tone effects 
because of the transformation of the multi-level picture 
into spatially encoded representations compatible with 
binary displays. A continuous tone image is one which 
contains an apparent continuum of gray levels. Some scenes 
wh en viewed by the human eye may require up to or more than 
256 discrete gray levels to give the appearance of a 
continuum of gray levels from one shade to another. 

Examples of continuous tone images are television (CRT) 

images, photographic images and real world scenes viewed by 
a camera. 

The bilevel/binary output images produced by the 
algorithms discussed in the subsequent chapters differ from 
other bilevel images like halftone images, and line copy 
images. In the bilevel pictures produced by these algorithms, 
the dots (or the dark spots) in the output have Ijhe same 
size and the gray level effect is obtained by their spatial 
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arrangement. Algorithms discussed for achieving this 
spatial arrangement are of two typesi one where the 
resolution of both the input as well as the output image 
is the same, and second, when the resolution of the output 
image is higher than the input image. 

Halftone imagery is also a bilevel image which gives 
the subjective effect of a continuum of gray levels. The 
halftone techniques were developed in the mid 19th century 
for approximating the gray levels available in the natural 
imagery. There are three main techniques - lithography, 
letter press and gavure? and all these produce the gray 
levels via the presence or absence of opaque ink on a page. 

In order to represent the continuum of gray levels, high 
frequency line and dot structures are printed which have their 
width varied spatially throughout the scene to yield a 
varying percent reflectance across the page. The end 
result is that when such images are viewed at normal 
viewing distances, the dot or line structure is not 
noticeable but the varying average gray level produces an 
approximation to a natural scene. Examples of halftone 
images are the pictures printed in books and newspapers, 

ITie. third class of bilevel pictures is the lihe copy 
imagery. This imagery is composed of alphanumeric ' 

characters^,, lib© - segwnts ,and solid areas of a 



single gray level. This class of images are essentially 
made up of only two gray levels but unlike the halftone 
images, the dots and lines of visible size are used. 

Examples of this imagery are books, magazines (less their 
halftone pictures) and maps, 

1*1 DISPLAY DEVICES 

Many displays are basically bilevel in nature which 
have individual display cells of the same size arranged in 
a rectangular array. Individual display cells or sites 
are either on or offj bright or dark> white or black) 
reflective or absorptive. Examples of these kinds of 
displays are * Direct view storage tube (DVST), liquid 
crystal displays, AC— sustained gas discharge or plasma 
panels. The algorithms studied are used for display of 
images on such displays. 

These display devices have inherent image storage 
capability and overcome the disadvantage of refresh displays 
of CRT, DVST resembles a CRT in that it uses a similar 
electron gun and a phosphorus coated screen. The difference 
is that the phosphorus has extremely long persistence and 
the beam does not write directly on the phosphoius,but on a 
fine mesh wire-grid coated with dielectric and mounted 
just behind the Screen. A pattern of positive charge is 
deposited on the grid, and this pattern is transferred to 
the phosphonashy a cdntinuous flood of electrons issuing 
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from a separate flood gun. Just behind the storage mesh 
is a second grid, the collector, whose main aim is to smooth 
out the flow of flood electrons. These electrons pass 
through the collector at a low velocity, and are attracted 
to the positively charged portions of the storage mesh but 
repelled by the rest. Electrons not repelled by the 
storage mesh pass right through it and strike the phosphorus. 
To increase the energy of these slow moving electrons to 
create a bright picture, the screen is maintained at a high 
positive potential. Fig, 1,1 shows the general arrangement 
of the DVST. The outputs of the study have been generated on 
such a device, the Tektronix 4006-1, which incorporates a 
7 by 10 inch DVST and a built in alphanumeric keyboard. 

The AC-plasma panel functionally is very similar to 
the DVST, even though its construction is very different. 

Images can be written onto the display surface point by 
pointj each point remains bright after it has been intensi- 
fied. Construction of the plasma panel is shown in Fig, 1,2, 
It consists of two sheets of glass with thin, closely 
spaced gold electrodes attached to the inner faces and 
covered with a dielectric material. The intervening space 
between the sheets of glass is filled with neon-based gas 
and sealed. By applying voltages between the electrodes the 
gas within the panel is made to behave as if it were divided 
into tiny cells# each #ne independent of its neighbours, A 
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cell is made to glow by placing a firing voltage across it 
by means of electrodes. The gas within the cell begins to 
discharge, and this develops rapidly into a glow. The glow 
is sustained by maintaining a high frequency AC voltage 
across the cell. 

Among the newer display techniques, the thermal effects 
in liquid crystals are made use of to display pictures. 

Fig. 1,3 shows an arrangement using liquid crystals. It 
employs a laser beam for writing. Here, all the liquid- 
crystal molecules are assumed to be initially aligned 
perpendicular to the glass plates, making the layer appear 
clear. When scanned with the modulated laser beam the 
energy absorbed results in local heating. This raises the 
material above the transition temperature, causing the 
molecules to become completely disordered. After the removal 
of the laser beam, the material cools rapidly and the 
molecules remain randomly oriented, producing a scattering 
effect which stays for very long. The antireflection 
layer is provided to increase the absorption of energy from 
the laser beam. Light which enters clear elements of the 
liquid crystal is reflected back from the rear aluminium 
electrode. This light reaches the screen through the 
aperture. However, the light reflected from the scattered 
elements is emitted in a variety of directions with a large 
fraction being unable to pass through the aperture. Since 
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Fig. 1.3 The liquid crystal display 
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the lens also serves to image the liquid crystal layer into 
the viewing screen, elements addressed by the laser ‘beam 
appear as a pattern of dark spots against a light background 
on the screen. The above devices are discussed in detail in 
[l] and [2]. 

1.2 SCOPE OF THE STUDY 

We can model the picture processing problem as given 
in Fig. 1.4. The input is restricted only to multilevel 
computer generated and continuous tone pictures and the 
output is a binary picture. The algorithms are used to 
process the pictorial input for high fidelity binary output 
pictures. The processing involves the selection of a 
threshold value which is used to spatially arrange the dots/ 
blanks at the output. 

The study has been carried out with the help of a 32- 
gray level computer generated picture of Abraham Lincoln. 

This picture has been chosen from a set of 32-gray level 
pictures given as 64x64 matrix of sampled intensity values 
at Appendix *A* of [3]. The algorithms- . .o-?- are applicable 
to any continuous tone picture whose scanned intensities are 
available. The processing has been done by simulation on the 
DEC- 10 system and the hard copies of Lincoln’s pictures have 
been obtained with the help of the plotter available with the 
Tektronix 4001-6 graphics terminal. The pictures generated 
are of three sizes*, 64x64, 128x128 and 256x256, 
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The available 64x64, 32-lGvel data of Lincoln’s 
picture was processed to generate the 128x128 and 
256x256 sized data. The data was processed in two ways : 
repetitive and random. In the former, each pixel’s gray 
value was replaced by a 2x2 or 4x4 matrix consisting of the 
same pixel gray value repeated to generate 128x128 and 
256x256 sized data respectively. For random processing the 
data was generated by generating a set of random numbers 
lying in the interval between the pixel's gray value and 
the gray average over a 3x3 window. The required number of 
these random numbers were then chosen and used to replace 
the original value. The computer programs used for 
generating these data, DIST. FOR and FOIST .FOR are attached 
at Appendix B. To generate the random numbers the NAG 
library subroutine G05DDF was used. 

To evaluate and compare the various algorithms, the 
correct method perhaps would be to assess the pictorial 
outputs by subjective tests involving a number of subjects 
evaluating and responding to the quality of the picture. 
Since that is a long process and will probably be more 
suitable for a specific application, three image quality 
metrics given below were chosen as reasonable components 
or embodiments of the subjective judgements. 

( a) Low Frequency Rendition : Capability of the process 
to reproduce the low frequency pictorial information in the 
input. 
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(b) High frequency Rendition s Capability of the process 
to reproduce fine details. That is to say that an algorithm 
capable of reproducing the. same contrast as that of the 
input will be superior to an algorithm which reproduces the 
same fine detail with lesser or no contrast, 

(c) Processing Artifacts j This implies those details in 
the output picture which are not part of the original 
picture but have appeared as a result of the processing 
technique. False contours are a type of artifact and may 
be caused due to the gray level quantization step which are 
sufficiently large to create a visible contour when the 
input image is truly a smooth, gradual variation from 

one to the other gray level. Another type of artifact 
which is similar to the false contour is caused by the 
output patterns. In order to get a gray level with the 
help of the binary process, called orthographic technique, 
use of binary patterns over some area is made to provide 
an average reflectance which is equivalent to a desired 
gray level. Such output patterns may create certain 
textures. 

Another metric which may bo used to compare and 
evaluate the various algorithms is the complexity of each 
of these algorithms in implementation. Though the algorithms 
are not very complex but keeping in view the real time 
implementation, even small computational difference may also 
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turn out to be of reasonable importance. The measure 
utilized is ttie ‘context’ required by each algorithm. 

The context is defined as the size of the memory required to 
determine the output state of a given pixel. For example, 
in a single level thresholding (like the fixed thresholding) 
the decision over the state corresponding to a given pixel 
is simply made by just that pixel being made available 
whereas, in case of an algorithm where a window (local) 
average has to be computed before the state of the output 
pixel can be decided, the context is the number of pixels 
required in the computation of the average. As seen, it 
does not consider operations like number of multiplications 
etc. to gauge the complexity. 

No study of image processing is complete without 
reference to the capability of the human visual system as 
all the pictures are ultimately viewed by humans only. Towards 
this end the perception by the eye is of specific importance 
as the eye sees the gray level in the binary picture due to 
its inherent property of spatial integration. Efforts have 
been made to approximate the eye by a mathematical model. 

One such model has been chosen from the various models 
available in the literature and the processing has been done 
for a set, of 3x3 binary patterns as inputs. The output 
2-dimensional arrays give the mathematical approximations to 
what the eye sees. Some of these outputs have been plotted 
with the help of 2-dimensional hidden line plots to give the 



12 


effect of a surface seen by the eye. This aspect has 
been covered in Chapter II. Subsequently sets of 4x4 
and 6x6 binary patterns have been chosen to output the 
Lincoln's picture by the orthographic technique. The 
resulting picture is a 256x256 and 384x384 pixel picture 
respectively as for each pixel of the input picture we 
print a 4x4 or 6x6 binary pattern at the output. 

Other techniques discussed differ from the 
orthographic technique in that they have the same resolution 
at the input and output. Eye’s perception of gray levels 
in such pictures is spread over a small area to give a 
general effect of an average gray level. 

To facilitate comparison between the various 

algorithms, the following notations have been used, A frame 

2 

is defined as a set of N samples taken on a square geometric 

grid of N lines with N samples per line. In this study, 

three different values of N (64,128 and 256)have been con- 

2 

sidered. Each of the N samples is the intensity value 
I(x,y) of the point (x,y) in the picture, subscripts x and y 
fixing the position within a line and within a frame. Since 
the Lincoln's picture is a 32 level picture, I(x,y) =0 
implies the darkest gray level and I(x,y) = 31 means the 
brightest gray level, 0(x,y) denotes the displayed inten- 
sity which is either 0 or 31 as the picture displayed either 
has a dot or a blank. 
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CHAPTER II 

THE EVE AND ITS MODELS 

Light as per Webster* s dictionary is : 'radiant energy 
which by its action on the organs of vision, enables them 
to perform their function of vision'. Light as we know is 
a form of electromagnetic radiation lying in the relatively 
narrow region of the electromagnetic spectrum over a wave- 
length band of about 350 to 780 nm. The mechanisms by 
which the light interacts with the 'organs of vision* is 
not very well understood. However, an attempt is made in 
this chapter to review some of the mathematical models of 
the eye available in the literature and then one of them 
is chosen and used to show analytically how the eye 
perceives a fine binary pattern structure(of dots and blanks) 
as different gray levels. 

2.1 THE EYE 

A conceptual technique for the establishment of a 
model of the human visual system (HVS) would be to perform 
a physiological analysis of the eye, the wave paths to the 
brain, and those parts of the brain involved in visual 
perception. Such a task, of course, is beyond the ability 
of the man because of the large number of infinitesimally 
small elements in the visual chain and their complex inter- 
connections. Fig. 2.1 contains a sketch of the horizontal 
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cross-section of the human eyeball. The front of the eye 
is covered by a transparent surface called the cornea. The 
remaining outer cover called the scalera, is composed 
of a fibrous coat that surrounds the choroid, a layer 
containing blood capillaries. Inside the choroid is the 
retina, which is composed of two type of receptors - rods 
and cones. Nerves connecting to the retina leave the 
eyeball through the optic nerve bundle. Light entering the 
cornea is focussed to the retinal surface by a lens that 
changes shape under muscular control to perform proper 
focussing of near and distant objects. An iris acts as 
a diaphragm to control the amount of light entering the eye. 
The rods in the retina are long slender receptors, while 
the cones are generally shorter and thicker in structure. 
Operationally, the rods are more sensitive than cones to 
light. At low levels of illumination, the rods provide a 
visual response called scotopic vision. Cones respond to 
higher levels of illuminationj their response is called 
photopic vision. An eye contains about 6,5 million cones 
and 100 million rods distributed over the retina. Fig, 2,2 
shows the distribution of rods and cones over a horizontal 
line on the retina. At a point near the optic nerve called 
the fovea, the density of cones is greatest. This is the 
region of sharpest photopic vision. There are neither rods 
nor cones in the vicinity of the optic nerve, and hence the 
eye has a 'blind spot’ in this region. When a light 
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stimulus activates a rod or cone, a photochemical 
transition occurs producing a nerve impulse. The manner in 
which the nerve impulses propagate through the visual 
system is presently not well established. It is known that 
the optic nerve bundle contains on the order of 800,000 
nerve fibres. Since there are 100,000,000 (lOO million) 
receptors in the retina, it is obvious that in many 
regions of the retina, the rods and cones must be inter- 
connected to nerve fibres. Since, neither the photo- 
chemistry nor the propagation of nerve impulses within the 
eye is well understood, a deterministic characterization 
of the visual process is not available. However, certain 
models of the HVS have been presented which are taken to 
predict the human visual system response. 

Before the HVS model is discussed, certain visual 
phenomena that a reasonable model must accomodate are 
discussed s 

a) Dark adaptation of the eye : 

It is well known that a person is more sensitive to dim 
flashes of light if he has been in dark for a period of time 
(as happens when one enters a movie hall). Cornsweet [l] 
has shown that the cones are adapted to dark in about 
5 minutes whereas the rods take much longer (about 30 minutes). 
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b) Threshold intensity for perception : 

It is the lowest possible intensity that the eye can see. 
It is dependant upon the colour, location with respect to 
the eye, timing and size of the light. The state of the 
subjects eye will also affect the threshold value. Any 
light lost before it reaches the rods (being the primary 
receptors at low intensities), cannot affect the subject as 
he does not assimilate any information about the light if 
it is not absorbed by the visual pigment in the rods. It 
has been shown by Cornsweet (p.25) that after catering for 
all the losses of light like reflection from the eye, absorp- 
tion by the pigment that fill the eyeball itself, and the 
loss due to the inter rod spacings; about 10 quanta are 
absorbed, but they are spaced over an area which has 500 rods. 
Probability of one rod absorbing two quanta of light will 
be very small and hence one quantum is sufficient to stimulate 
a single rod and approximately lO rods within a 10 pm diameter 
area must be stimulated within 1 msec of each other before 
the stimulus is perceived. Hence, the activation of 10 rods 
is somehow added up in the visual system. This threshold 
is for monocular vision. The binocular threshold is lower 
by a factor of "{2 as mentioned by Hall [2], 

c) Simultaneous contrasl^ 

It is commonly experienced that the brightness of a 
region does not depend upon its intensity only but also upon 
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the surround. If we look at a thick (opaque) gray or 
coloured paper with one eye under any ordinary light and 
then hold it up at an aim*s length between the eye and 
a bright sky or a large light source, the paper will look 
much darker when it is in front of the light source, 
even though the actual amount of light reflected from the 
paper is the same. In fact, if the background light is 
bright enough the paper will look black when held in front 
of it. 

d) Brightness constancy s 

It implies that the brightness of any object remains 
fairly constant despite large changes in the illumination 
that falls on it. For example, a piece of ordinary white 
paper reflects about 99 percent of the incident light, while 
black papeis reflects about 10 percent, When these papers 
are moved from ordinary room light illumination into direct 
sunlight, the incident illumination and therefore, the 
reflected light intensity may increase by 1000 foldj yet 
the papers still have about the same brightness as they 
had indoors. Furthermore, the black paper outside reflects 
0.10x1000 s= 100 relative units of intensity and the white 
papers inside reflects only 0.9x1 = 0,9 units* but the black 
paper outside still looks darker than the white paper inside. 
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e) Mach bands : 

Consider a set of gray scale strips, each having the 
same gray level, the darkest strip towards the left extreme 
and the brightest towards the right extreme. The reflected 
light from each strip is uniform over its width and differs 
from its neighbours by a constant amount. Nevertheless, the 
visual appearance is that each strip is darker at its 
right side than at its left side. This is called the Mach 
Band Effect. The effect is due to the spatial frequency 
response of the eye, which possesses a lower sensitivity 
to high and low frequencies than mid-frequencies* 

2,2 MODULATION TRANSFER FUNCTION (MTF) 

In order to postulate a model for HVS, it is necessary 
to have a method \Miich can be used to predict the response 
of the HVS given a particular input. Since our concern 
in this study is the spatial response of the eye, the 
processing model is based upon the amenable frequency 
measurements. 

Since any spatial pattern can be broken down into 
various sine wave components by Fourier Analysis, the 
procedure to measure the MTF is to feed a sine wave grating 
as the input. Cornsweet [l] has explained this method with 
the help of an optical lens to which sine wave gratings of 
different frequencies are input and the outputs are studied 
to calculate the MTF of the lens. It is seen that as the 



21 


frequency of the gratings is increased, the amplitude of 
the modulation of the image distribution (i.e«, the 
difference between the intensity at each bright part 
and at each dark part) becomes smaller and smaller and 
for extremely fine gratings, the modulation in the image 
becomes zero and a uniform intensity distribution is seen. 

The property of the lens being measured is its ability to 
transfer spatial modulation of intensity from the object to 
the image. Hence, the MTF of the lens is an indication of 
its ability to resolve gratings. In the experiment above, 
the phase shifts have been assumed to be zero. If say, 
the phase shift introduced by the lens was constant for all 
the frequencies, then the image would merely be shifted 
from the expected position. However, if the phase shifts 
vary with frequency then the MTF alone would not be sufficient 
to define the image of any object as the relative 
positions of the various sine wave components would not be 
known. Therefore, if there were substantial phase shifts 
varying with frequency, both the MTF and the phase shift at 
each frequency would have to be known in order to predict the 
image correctly. However, the phase shifts produced by most 
lenses and by the optical system of the human eye are relati- 
vely small near the optic axis (fovea for the eye) and only 
trivial errors are introduced by assuming zero phase shifts. 
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Before we use the concept of MTF for the eye model, 
it is necessary and essential to list out the conditions 
that must be met in order to use Fourier techniques* 

a) Linearity : 

A linear system is one that obeys the principle of 
superposition, i*e,, if two inputs give the outputs X and Y, 
then when both the inputs are applied simultaneously, 
the system output should be X+Y. In general, if the 
intensity radiated from the object is increased, the 
magnitude of the response of the system should increase 
proportionately* We may also apply the technique over the 
linear operating region of a nonlinear system. 

b) Homogenity j 

A system is spatially homogeneous if its characteristics 
are the same in all locations* If a system is homogeneous, 
then shifts in the location of the input pattern may cause 
shifts in the location of the output pattern, but the 
output pattern will not change except in position. The HVS 
is not homogeneous due to the variation of the densities of 
rods and cones with retinal position. However, it is 
relatively homogeneous near the optic axis (i.e, near the 
fovea). In spite of the inhoraogenity of the HVS, Fourier 
techniques may be used to correctly predict its response. 
Cornsweet [l] has pointed out that there is probably a self 
homogenising process in the structurally anatomically 
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inhomogeneous process which renders it open for the use 
of Fourier techniques. Models discussed do not consider 
temporal responses, hence temporal homogenity has not 
been considered. 

c) Isotropic: 

A system is isotropic if its characteristics are 
the same in all directions. For example, a system is 
isotropic if the MTF measured with sinusoidal grating 
objects is the same regardless of whether the gratings 
are vertical, horizontal or oblique, (i.e, regardless of 
their angular orientation with respect to the optic axis 
in the object plane). The response of the HVS to a rotated 
contrast grating is a function of frequency of the grating 
as well as the angle of orientation. The sensitivity of the 
system decreases to a minimum at 45° and then rises again 
reaching the original level at 90°, At the point of maximum 
deviation (i.e, 45°) a frequency of 30 cycles/degree would 
be -3 dB below the response at 0° rotation. Whereas for 
10 cycles/degree the response at 45° deviation is only 
15 percent lower than the original. This an-isotropic 
behaviour has not been included in the models discussed. 

2.3 THE EYE MODEL 

In the lens experiment explained in Section 2.2, the 
output to a particular input is visible and the transfer 
function can be calculated in a straight-forward manner. The 
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case is different for the eye, where we have no way of 
knowing what kind of output image is being formed in 
the human brain# Hence, in order to measure the modulation 
transfer function for the HVS, various ways have been used. 
But the basic principle used in all is that a subject is 
asked to match a variable object (grating) with a reference 
object (grating) and this can be done in different ways 
like varying the contrast or frequency. 

As a direct analogy with the procedure for measuring 
the MTF of a lens, a human subject could be presented with a 
series of sine wave gratings with same amplitudes but 
different frequencies and then obtaining a measure of his 
perception of the patterns. The measure of the perception 
is the apparent amplitude of each sine wave, (that is, the 
difference between the brightness of peaks and troughs of 
the wave ) because the apparent amplitude is exactly an 
'output* of the HVS. Thus, a plot of the input frequency 
( for a fixed amplitude) against the apparent amplitude, 
would be analogous to the plot of the MTF of the lens. 

Given such a plot and assuming that the HVS is operating 
in the linear region, and that it is homogeneous and 
isotropic, wo should be able to predict the appearance of 
any input pattern, Davidson [l, p. 335] has suggested 
one such arrangement. Assuming the HVS to be linear, 
homogeneous, isotropic, monocular, monochromatic and 
photopie (that is operating beyond threshold) , a MTF was 
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plotted and is given in Fig, 2.3. In fact, the response 
is contrast dependant and the curves for different ratios 
of contrast have same shape with higher ratio curves, 
higher in amplitude. This response can be represented by 
a system which is just a band pass filter but it would be 
difficult to defend from a physiological standpoint since 
the response is a compound one due to several mechanisms 
within the HVS, 

A slightly more detailed model would be a combination 
of a LPF and a HPF. In fact, the complete response of the 
HVS as in Fig. 2,3 can be depicted by a variable sine wave 
grating of Fig. 2,4, The figure shows sine wave gratings 
with frequency increasing from left to right and the 
contrast increasing from bottom to top. If a line is 
drawn horizontally across the figure, the line will intersect 
sinusoidal changes in intensity whose amplitude is constant 
but whose frequency varies. As can be seen clearly from the 
diagram, sensitivity of the HVS is maximum at raid frequencies 
and reduces at lower and higher spatial frequencies, 

Cornsweet (p, 343) has used this MTF to explain the Mach band 
effect and contrast sensitivity mentioned earlier in this 
chapter. 

Characteristics of the high frequency response, i,e., 
the region that extends from peak mid frequency 

where the sensitivity decreases with frequency. 
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are determined in part by the optics of the eye (lens» 
pupil size and cornea)? also the properties of the retina 
itself affect the high frequency response. Size and density 
of the photo receptors, neural summation and scattering of 
light at the retinal surface affect the high frequency 
response of the HVS. The low frequency response is limited 
by the lateral inhibition. 

This model of a combination of LPF and HPF can be used 
to approximate a HVS but it does not account for the pheno- 
menon of ’brightness constancy'. This phenomenon can be 
accounted for by introducing a nonlinearity in the system. 
This is illustrated in Fig,, 2, 5(a) where the intensities 
in the right hand figure are all the same multiple of those 
on the left. If a logarithmic transformation is assumed, 
the figures get transformed as in Fig. 2.5(b). The two 
outputs now look quite similar which is how the eye perceives 
the black and white papers in normal illumination and under 
bright illumination source. That is, if a logarithmic trans 
formation is added to the MTF described earlier, we are in 
a position to accomodate in the model the majority of the 
phenomena of the eye's perception of still, monochrome images 
This nonlinearity has been studied extensively by various 
researchers. The studies were mainly carried out on the 
horseshoe crab, 'Limulus', and with the use of mechanical 
receptors. The choice of Limulus was motivated because 














the lateral eyes of this animal lend themselves to 
physiological analysis, and preliminary studies suggested 
that many of the properties of this eye are similar to 
those of higher animals, including humans. It has been 
shown that the nonlinearity in the HVS can indeed be 
approximated by a logarithmic function. 

The question which arises is how a nonlinearity 
can be added to MTF, (a linear system). It has been 
argued by Cornsweet (p, 334), that for small variations in 
intensity of input patterns, the system essentially operates 
in the linear region of the nonlinearity* Also, it has been 
shown that under increasing light adaptation conditions the 
linear region becomes larger. 

Two Ways in which the nonlinearity can be added are 
shown in Fig, 2,6. Of the two structures shown Fig. 2, 6(a) 
can be used to explain the phenomenon of brightness con- 
stancy, This system is thus physiologically sound and it 
also predicts the brightness constancy. This model has been 
chosen to process the 3x3 binary patterns. The impulse 
response of the LPF has been given by Campbell [3] as s 

h^(x) = exp(-a jx |) (2.1) 

where a has a value of 0,7 for a pupil diameter of 3 mm. 

Its Fourier transform is given by 
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The -3 dB point occurs at approximately 6*5 cycles/degree. 
The transfer function os HPF is given as j 


Ho(«) 


2 

a + 


w 


2aoa 


(l-aQ)(a^+ 


(2.3) 


where is the distance factor between receptors and *a' 
is the inhibiting factor. This transfer function has been 
derived by using the backward inhibition model for the 
photoreceptors in HVS given in Fig. 2.7. The receptors are 
assumed to have a logarithmic response to the incoming 
light. The principle used is that there is an inhibiting 
effect in between the receptors in the retina. The detailed 
analysis has been done in the work by Hall and Hall [2]. 

2.4 RESPONSE TO BINARY PATTERNS 


After having discussed a model for the HVS, the next 
step is to use it to analyse its response to a set of patterns 
of binary pixels. As was mentioned briefly in Chapter I, due 
to the spatial integration properties of the eye, a set of 
properly chosen closely spaced binary patterns with their 
basic repetition frequency beyond the eye’s cutoff range are 
perceived as different levels of brightness. This property 
is usefully utilized to transform multilevel computer genera- 
ted (or continuous tone) pictures into binary pictures which 
retain the gray level effects due to the spatial arrangement 
of the binary patterns. 
















32 


A set of 3x3 such binary patterns which can be used 
to display 10 intensities are given in Fig. 2.8. This 
set of binary patterns has been taken from the book by 
Newman and Sproull [4]. These patterns are input to the 
model of the HVS and the outputs obtained with the help \ 
of the computer program PROC«FOR. 

To find out the MTF, the impulse response of the LPF 
is convolved with the response of the HPF to give the overall 
impulse response of the model less the nonlinearity. The 
impulse response of the HPF is worked out as j 

^2^^^ ” i? ®xp(~D|x|) - ^ [D^ exp(-D)xj) - 2D 6(x)] (2.4) 

0.8 

0.005 

0.000125 

D (B+a^)^= 0.0714 
On substituting these values, we get 

hgCx) * 1.25 6(x) - 0.04375 exp(-0.07l4 Ixj) 

Hence, the overall response of the system is 

h(x) ss h^(x) » h2(x) = 1.2531 exp(-0.7 j x| ) 

-0.06125 exp(-0.07l4 (xf) 


(2.5) 


( 2 . 6 ) 
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Assuming circular symmetry, the two-dimensional response 
of the system is 

h(x) = 1.2531 exp(-0.7 jx+y 1) - 0.06125 exp(-0.0714 ix+ y I) 

(2.7) 

The input to the eye model is one of the 3x3 binary 
patterns in which the dot is assigned value 2 and blank is 
assigned value ll. The pattern has been transformed into 
a tv/o-dimensional array with the centre pixel being the 
origin.- In order to find the response of the system, we 
consider only the discrete values of the system point 
Spread function at the points at which the patterns exist. 

The systems response has been limited to a size of 2lx21 as 
beyond this region the values of the output are quite small. 
The system response thus is the two dimensional convolution 
of the log of the 3x3 pattern and the MTF. The two dimen- 
sional convolution has been tried out in two ways s the 
usual double summation and a method suggested by MacAdam [5] 
to convert the 2-dimensional convolution into 1-dimensional 
convolution. In this method the 2-dimensional arrays are 
read as vectors by the computer and then 1-dimensional 
convolution is performed to yield a vector which is rewritten 
as a 2-dimGnsional array to give the result. Though while 
using the DEC- 10 system the difference in the CPU time was 
negligible as the arrays were not large. 
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The computer program used to carry out this operation 

is named CONVOL*FOR and is attached at Appendix B. The 

program generates the 21x21 matrix which is the MTF of the 

system. The input to the program is the 3x3 binary pattern. 

The output is a 23x23 two dimensional array which has the 

centre value i.e., the 12th column of the 12th row 

maximum. These values progressively decrease and arc 

minimum at the four corners of the array. As seen from Appx.C for t 

10 patterns (and their combinations) the centre value is 

highest for the brightest level, i.e., the pattern having 

all bright squares and is lowest for the pattern having all 

dots, thus the 10 gray levels are distinguished as different 

from each other. To get a better idea of these outputs, 

some of these have been plotted with the help of 3-dimensional 

plots with and without hidden line elimination and are shown 
as Figs. 2.9 and 2.10. 

The 3-dimensional plotting without hidden line elimina- 
tion is done by fitting a curve (using cubic spline fit) 
for the values generated by CQMV0L*F0R. The algorithm used 
for this purpose has been suggested by Cline [6], Cline has 
suggested curve fitting by cubic splines with tension. The 
tension factor removes the defect of inflection points 
commonly seen in the cubic spline fitted curves. With the 
tension factori the curve may be assumed to be a light and 
flexible bar not only constrained to pass through the given 
values but also able to respond to a tension produced by 
pulling on its ends. Sufficient tension would remove the 
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unwanted inflection points. The program consists of 
the driver program SPLINE -FOR and two subroutines KURVl and 
KURV2. These are attached at Appendix B. The plot drawn is 
shown at Fig. 2,9, This does not incorporate the hidden 
line elimination. The curves have been rotated by 20° to 
give the 3-dimensional effect. 

The plotting with hidden line elimination is carried out 
with the help of an algorithm suggested by Watkins [7‘]. 

This algorithm incorporates hidden line elimination and the 
plots are rotated to give the desired 3-dimensional effect. 
The input to this program is the output of the SPLINE ‘FOR 
program. Program PLOT 3D- FOR accepts the three dimensional 
data, rotates it in the three space and plots the projection 
of the resulting figure onto the X-Y plane. One point to 
note here is that the rotation carried out by PLOTSD-FOR 
is not in addition to the rotation carried out by SPLINE- FOR 
because the data output by SPLINE* FOR is in unrotated form 
and the rotation that was achieved was with the help of the 
GPGS package while drawing the plot on the screen. 
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CHAPTER III 

TECHNIQUES FOR DISPLAY OF MULTILEVEL GRAY 
PICTURES ON BINARY DISPLAYS 

3.1 INTRODUCTION 

As mentioned in Chapter I, there are two classes of 
techniques to display multi-gray level (continuous tone 
or computer generated) pictures on binary displays. They 
differ in the relationship between the number of pixels 
of the input and output pictures. Gray levels were 
approximated by patterns of binary pixels in the previous 
chapter and the technique using this principle, called the 
Orthographic Technique, has higher number of pixels in the 
output picture as compared to the pixels in the input 
picture. The remaining techniques discussed in this chapter 
belong to the general class of algorithms that convey gray 
level information by the spatial arrangement of binary 
pixels which are in a one to one correspondence with the 
pixels in the original image. These techniques create 
patterns of binary pixels that preserve detail and average 
image brightness over extended areas, 

3.2 ORTHOGRAPHIC TECHNIQUE 

This technique utilizes an 'mxn' array of binary 
pixels in the form of a gray scale ‘character* to represent 
a multi-gray level picture. These characters together form a 
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gray scale *font’ which, when printed with minimal inter- 
character spacing, yield the gray scale infoimation of the 
original picture. 

In Chapter II a set of 3x3 binary patterns were shown 
and they could be used to represent a 10 gray level picture. 
Similarly, Hamill [ij has suggested a set of 4x4 binary 
pixel patterns depicting 16 gray levels. This set of 
patterns has been used to generate the test picture of 
Lincoln shown in Fig, l of Appendix A, The generated 
picture is shown at Fig,2a of Appendix A. The program 
used for the technique is called ORTHO- FOR and is attached 
at Appendix B* The algorithm simply involves a table look 
up process as shown in Fig. 3.1. Since the patterns 
number only 16, a quantization step of two has been used 
to represent the 32 gray levels. The program ORTHO* FOR 
reads the patterns as vectorsof 16 and then arranges 
them as 4x4 patterns at the output for each input pixel. 

The picture generated is thus a 256x256 pixel picture. 

Since the usual pictures are sampled with the help of 
128 , 256 or even higher gray levels, a compromise is usually 
developed for the size of the patterns as large sized patterns 
produce coarse resolution. In view of this, large quantiza- 
tion steps will be needed which result in tone scale errors. 
Because of the coarseness of tone scales, false contouring 
might result and the actual characters used to create gray 
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Fig. 3*1 Signal flow diagram for Orthographic 
Technique 
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Fig, 3.2 The basic processing organisation 
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scales may also become visible. The 4x4 patterns used 
for the generation of the picture in Fig. 2a of Appendix A 
are visible. In certain cases, depending upon the selection 
of fonts, texture contours might be visible in the output 
picture. There is no processing complexity as the technique 
simply requires a pixel look up table. 

The set of patterns suggested by Hamill were modified 
into 6x6 binary patterns to represent 32 gray levels. The 
picture generated by this set is as shov/n in Fig. 2b of 
Appendix A. Using the same set of 6x6 patterns another 
picture was generated and is shown at Fig, 9 of Appendix A. 
The quality of this picture is poorer than the picture of 
Fig, 2b of Appendix A. This was caused due to a programming 
error. A similar error in the program using Hamill' s 4x4 
patterns resulted in the absolutely distorted picture of 
Fig, 8 of Appendix A. In this, the data has been compressed 
and repeated to give the distortion where smaller sized 
pictures of Lincoln have been superimposed on the original 
picture, 

3,3 TECHNIQUES WITH EQUAL INPUT AND OUTPUT RESOLUTION 

The basic processing organisation that produces the 
binary output for the generation of the output picture is 
shown in Fig, 3,2. Intensity of each pixel of the input 
picture is compared with a threshold signal. If the 
intensity of the pixel is greater than the threshold the 
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corresponding output pixel is set to bright state, other- 
wise the output pixel is set to dark state. The various 
techniques described in this section differ primarily in. 
the way they produce the threshold values for comparison 
with the intensity of the input pixel. 

3.3,1 Fixed Threshold Technique 

Fig, 3.3 shows a schematic diagram for a globally 
fixed level thresholding process with input I(x,y) and 
output 0(x,y). This technique is simplest to implement 
in that the decision rule is simply stated as 

If I(x,y) > T; then 0(x,y) = 1 

else 0(x,y) = 0 

The threshold T is predecided and globally fixed and 
is generally in the neighbourhood of the centre value in 
the tone scale (16 for the test picture of Lincoln). The 
output process simply generates black and white pixels, 
depending upon the gray level of the input image at that 
location. Fig, 3a of Appendix A shows a picture generated 
by this method. As is obvious, the two ends of the gray 
scale are reproduced accurately but the remainder are 
severely distorted. The output picture either has dark or 
bright areas. 

Fine detail which has gray level swings above and 
below the threshold will be seen. Hence, mid range fine 
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o(x,y) 


Fig* 3.3 Signal flow diagram for globally 
fixed threshold 
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Fig. 3.4 Dither matricos (a) proposed by Lippel and 
Kurland (b) proposed by Bayer 
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details is reproduced but the details away from the mid- 
range are never reproduced. Though the detail rendition 
is limited to a very small range, very high frequency 
details are enabled in that range, A high contrast text 
can be displayed effectively by this technique but not a 
multi-gray level picture. The processing complexity of 
this technique is minimum and only one pixel of 'context’ 
is required for the computation. 

Fig. 3b of Appendix A shows a picture that has been 
generated in a slightly different manner. Instead of 
comparing a pixel of input picture, an average of the pixel 
and its eight immediate neighbours is calculated and 
compared to the fixed threshold. As seen, the quality of 
the output picture still remains the same. The processing 
complexity of this technique is more as nine pixels of 
’context' are required to compute the window average. 

3.3.2 Ordered Dither Technique 

This technique was worked out for TV signals which 
being analog had to bo quantized before digital transmission 
and due to the need to reduce the number of bits trans- 
mitted, the obvious solution was to reduce the number of 
quantization steps by increasing the step size. This 
resulted in the degradation of the TV picture in the areas 
where intensity changed slowly (i.e. low detail areas). The 
degradation took the form of 'false contours’. To eliminate 
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these contours it was suggested to add a pseudo noise 
signal called a dither signal to the true input signal* 

This produced rapid switching between the quantizer levels 
on either side of the true input signal* Pictures were 
transmitted by using as few as 4 bits/pixel (i*e*, 16 gray 
levels) by Roberts [2] by adding a 2 dimensional pseudo 
noise sequence* The false contouring in the received 
picture was reduced by the fact that the subject viewed 
the picture as a near continuum of gray levels by the 
eye's averaging over a small two dimensional area. 

In the present context the multilevel picture is being 
quantized into just two levels and one may add a two 
dimensional pseudo noise sequence to an input prior to 
quantization to two levels. The resulting picture has the 
errors distributed spatially and the observer integrates 

i 

the average reflectance in a small region and sees a near 
continuum of gray levels. Limb [3], Lippel and Kurland [4] 
have used it for binary representation* It has been shown 
subjectively that 'ordered dither patterns' constructed with 
a small square dither matrix repeated over the picture area 
to form a rectangular array of matrices is superior to a 
random arrangement of the same population of samples* 

Pictures generated by the two kinds of dither have been 
studied by Lippel and Kurland [4], They have also studied 
the arrangement of dither samples in the dither matrix to 
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find out optiraum dither patterns keeping in mind the effect 
the pattern has on pictorial information in the output 
picture* The dither patterns suggested by them are given 
in Fig* 3.4(a)* Bayer has also suggested a dither pattern 
given in Fig* 3.4(b). 

This technique thus generates a bilevel representation 
of continuous tone images by comparing the image values 
I(X|,y) to a position dependant set of thresholds contained 
in a square dither matrix. This dither matrix is repeated 
over the entire picture in a checkerboard fashion and the 
decision for a pixel to be turned bright or dark is 


If I(x,y) > D(x,y), then 0(x,y) = 1 

else 0(x,y) = 0 


Jarvis et al [5] have shown a recursive relationship 

between the 2x2 matrix of Limb and the matrix generated by 

2 

Bayer* The Limb’s matrix D is one of the four 2x2 matrices 
that satisfy the Bayer optimization criteria* 
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This matrix scaled by a factor of 2 to cover the entire 
range of 32 levels has been used to generate the picture of 
Fig. 4a and Fig. 4c of Appendix A. Fig. 4a is a 128x128 
picture whereas Fig. 4c is a 256x256 picture. Bayer’s 
matrix of Fig, 3.4b has been used to generate the picture 
of Fig. 4b of Appendix A. The pictures of Figs.4d, 4e and 
4f of Appendix A have been generated using the matrix of 
(3.4) above but without any scaling to cover tho entire 
range of tho tone scale. 

This technique results in complete gray level scale 
generation in the output picture . Comparing pictures of 
Figs. 4a and 4c of Appendix A with the picture generated by 
Bayer’s matrix (Fig. 4b), it is seen that Bayer’s matrix 
produces more detail. It is due to the entire gray scale 
range being covered in steps of i as compared to the 
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steps of two for the pictures of Fig. 4a and 4c. The 
pictures of Figs. 4d, 4e and 4f miss out half the tonal 
scale due to error in the matrix used. However, picture of 
4f gives slightly more detail due to the randomness of the 
data used. 

The processing complexity of this technique is slightly 
more than ttie fixed threshold technique as here also only 
one pixel is processed at a time and instead of a globally 
fixed threshold it has a look up table for the thresholds, 

3,3#3 Constrained Average Technique 

This technique developed by Jarvis and Roberts [6] 
generates one bit for each pixel of input picture I(x,y) by 
comparing its intensity to a threshold value T computed from 
I(x,y) and its eight nearest neighbours in the array, of the 
input pixels. If I(x,y) is greater than the threshold T, 
the corresponding pixel at the output is turned bright. This 
technique provides edge emphasis alongwith gray scale rendi— 
tion* The process uses the finite signal to noise ratio 
inherent in the picture data due to sampling and quantization 
errors, to generate, on a statistical basis a number of 
bright cells related to the average image intensity in the 
picture area# 


The threshold is computed as : 
T « Y t I(x,y) [l - 


( 3 . 5 ) 
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where T{x,y) is the average intensity of I(x,y) and its 
eight nearest neighbours, R is the highest gray level 
(31 in the case of the test picture of Lincoln) and y is a 
parameter which controls the apparent contrast of the output 
image and is related to noise statistics. Since I(x,y) 
varies from 0 to R, the threshold changes in the value from 
y ro R-y, 

The working of the technique can be 
explained by the following argiments. 

When the average intensity in an area is equal to r/ 2, 
the threshold value from (3»5) is also equal to r/ 2 and an 
approximately equal number of pixels will have values above 
and bolow the threshold T resulting in one half of the 
pixels being turned bright* As the brightness increases, 
the threshold T is less than the average brightness T(x,y) 
resulting in more than one-half of the pixels being turned 
bright. The reverse happens when the brightness decreases, 

Tho techniques capability of enhancing edges is seen 
whon a pixel markedly different from the local average 
intensity is considered. The pixel will be turned on 
or off doponding upon its intensity value relative to the 
neighbouihood average. Pixels near a dark line (edge) will 
have higher probability of being turned on as the neighbour- 
hood average will be reduced due to the presence of the dark 
pixels Je«>It»«aAna to the line. The reverse will hold good 
for the pixels on "ttve dark lino itself as then the 
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neighboujrhood average would be more than the value of the 
pixel and it will be turned off. This suggests that the 
size of the area considered for the computation of the 
neighbourhood average and the threshold should be small 
as any step intensity changes or gradients in the area 
effect the state of the output pixels. Roberts [6] has 
shown that tho constant y depends upon the noise chara- 
cteristics of the signal. Thus, a usable value of y can 
bo obtained if we know the noise characteristics of the 
signal being processed. In practice, y may be thought of 
as an arbitrary parameter that can be specified to 
obtain a pleasing rendition of the original picture, y 
re suite In normal gray scale rendition. 

Fig. 3*5 shows a signal flow diagram for the con- 
strained average threshold technique. Figs. 5a to 5e and 
6a and 6b of Appendix A show the various pictures generated 
with different types of data, and different values of y. As 
can be seen, edge emphasis has been achieved in all the 
pictures irrespective of the value of y or the type of data 
(repetitive or random) used. For y = 6, the picture 
gonoratod only emphasizes the edges whereas pictures with 
Y > 0 render edge emphasis alongwith the gray scales. The 
control over the apparent contrast of the technique is 
visible in pictures 5b and 5c. Picture 5c has more contrast 
than 5b. With y < 0, the generated picture is a negative 

j.l.T. kantor F 
l ^rgAL ubraW 
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3,5 Sifuia.l processing flow di^aram for 
consLrainod Average technique 


E(x,y) 



Fig. 3.6a Signal flow diagram for the error 
diffusion method 
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Fig. 3.6(b) The size of the error matrix 
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like picturo but it is not a true negative as the edge 
enhanceiBent still turns on cells that are more positive 
than their neighbourhood average* Finally > picture of 5f 
has been generated by adding a dither signal having 
rectangular distribution to the data of the picture* This 
gives the bast effects of gray scales and maximum details. 
The complexity of the constrained Average Technique is more 
than both the fixed threshold and ordered dither techniques,i 
as a ‘context’ of 9 pixels is required to compute the 
neighbourhood average and the technique requires further 
processing to calculate the threshold. 

3# 3* 4 Error Diffusion Technique 

This technique generates local patterns of binary 
pixels in such a way as to minimize the intensity error 
between the source and displayed image. When a pixel which 
was initially gray is turned to bright or dark, it con- 
stitutes an error in gray scale. Thus to correct the 
picture, this gray error is redistributed to neaitoy pixels, 
restoring the total gray content to its original value. The 
original work was presented by Schroeder [7] where the 
pictures were processed with the help of more than two 
output states. The technique was used for binary displays 
by Floyd [8]. Given the gray level of input pixel, the 
binary decision about the output pixel is made and the error 
between these two levels is dispersed (diffused) to the 
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riobt snd bolow ■the processed pixel. Fig, 3,6a shows 

a signal flow diagram representation of the error 

diffusion technique. 


At a point (x,y) in the picture there will be an error 
between the displayed intensity and the original image 
which contains more intensity levels than the displayed image. 

E(x,y) = I(x,y) - 0(x,y) 

whore E(x,y) is the intensity error at a point (x,y), 

0(x,y) is the displayed intensity and I(x,y) is the original 
image intensity. To implement the algorithm, a corrected 
intensity I*(x,y) is computed from the previously computed 
errors and the intensity I(x,y) s 


I(x,y) = I(x,y) + I a.. E(i,j) 

ij 


In this equation, the ranges of the indices i and j are 
defined by the definition of the neighbourhood of point 
I(x,y), define the weights that multiply the gray 

scale errors E(i,j) for those pixels prior to the addition 


of gray scale for I(x,y), Fig. 3.6b gives the neighbourhood 
of the point (x,y) and the weights that multiply the 
various E(i,J). The corrected intensity is then compared to 
a fixed threshold (centre value of the tone scale) to deter- 
mine whether the corresponding display element should be 
bright or dark. That is, • 
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If I’(x,y) > Ti Then 0(x,y) = 1 

else 0(x,y) = 0, 

Finally the error is computed as 
E(K,y) = I*(x,y) ~ 0(x,y). 

The element marked * in Fig. 3.6b is the position of 
the point I(x,y) under consideration and does not enter 
into the error computation. Likewise, the elements marked - 
are ahead of thu point (x,y) in the scan line and no error 
for theso elements exists j hence, they can not enter into 
the error computation. Jarvis et al [5] have shown 
empirically ttiat the minimum size of the matrix (neighbourhood) 
needed to generate subjectively satisfactory results is of 
the order shown in Fig. 3.6b, Hence, the context of this 
technique is the order of (ij). In the present case, 

’context* is 12, as twelve error values are needed to 
compute the binary output. The complexity of the technique 
can also be visualised from the fact that it requires 3 line 
memories to contain the successive error values and also 
a substantial number of computations for each pixel are 
needed, as a weighted sum is required to be done. The 
boundary handling consists of setting the error values to 
zero for a boundary two pixels deep on each side and the 
top of the picture. Only the pixels inside this region are 
evaluated. This two pixel width boundary is normally not 
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vj.si.bX6 in q 128 x 128 ot higher picture* 

This technique, like the ordered dither technique, 
is capable of reproducing the complete tone scale but 
the spatial frequencies in the output image will be a 
function of both the gray level being reproduced and the 
input imago detail. This is due to the effect of the error 
diffusion to the neighbourhood. For example, the spatial 
distance between black pixels in a bright area will be 
very long as compared to that distance for mid- tone gray 
levels (due to the difference in the quantum of error in 
the two cases). There may also be a spatial translation in 
the location of an edge detail depending upon the image 
content above and to the left of the edge. 

This technique is capable of creating output signals 
which have lower spatial frequencies (as mentioned above) 
and arc thus more visible than in the case of ordered 
dither technique, which puts a limit on the low frequency 
response. 

Figs. 7a to 7c of Appendix A are the pictures generated 
by this technique. Picture of 7a is a 64x64 picture whereas 
the remaining pictures are 128x128. The picture of Fig. 7c 
has been generated from the random data and hence has more 
noise as compared to the picture of Fig. 7b which has been 
generated from the repetitive data. 
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CHAPTER IV 
CONCLUSION 

Electronic manipulation of the pictorial imagery 
for storagOf communication and reproduction often use 
binary marking/display processes for cost effectiveness. 
The binary data generated by these processes may be 
further coded to affect more saving. To illustrate, 
the picture of Lincoln after being converted to binary 
bits is Run Length Coded (RLC) and again reproduced. The 
program used to carry out this processing is named RLC. 

FOR and is attached at Appendix B. The picture of Lincoln 
generated by this program is attached as Fig. 10 of 
Appendix A. The picture is an exact replica of the 
constrained average picture of Fig. 5b of Appendix A 
which was taken as the binary picture for coding. The 
saving effected by the constrained average process was 
to reduce the data to l/5th of its original size and the 
run length coding fcirther affected a saving of slightly 
loss than half. 

Though the possibility of displaying pictures with 
full range of gray scales as a series of on or off cells 
on a grid of moderate resolution appears marginal, the 
pic ture* reproduced in this study show that usable pictures 
can be displayed in this fashion and hence provide a range 
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of options to tho designer of a visual communication system 
using bilevel displays. 

The processed pictures generated by the various 
algorithms have been explained in each of the sections/ 
sub— sections under which the algorithms have been discussed. 
The techniques are evaluated as under : 

4.1 EVALUATION OF TECHNIQUES 

(a) Orthographic Technique j 

The frequency rendition is restricted by the factors 
of m and n as the binary patterns used are of the size 
mxn. It is quite common to find false contouring caused due 
to the binary patterns used. It can be seen to a little 
oxtont in the picture of Fig. 2 of Appendix A. It will be 
more prominent when a higher gray level picture is reproduced 
by this techniouG. The output picture also contains a 
texture depending upon the patterns used. The binary pattern 
used can easily be seen in the picture of Fig. 2 of 
Appendix A. 

(b) Fixed Threshold Technique s 

The frequency rendition of this technique is restricted 
to only the mid-range gray levels. However, very high 
frequency fine details are rendered in this range. The 
black/white boundary in the output picture is an ever 
present and most common artefact (false contour) to be seen. 
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Low noise riding the gray levels near the threshold value 
gets aniplitiodo This low noise inay* in fact^ be due to the 
inherent noise characteristics (practical limits on signal to 
noise factor) of the scanning systems. It may result into a 
visible artefact which was not present in the original 
picture. 

(c) Ordered Dither i 

Due to the checkerboard kind of arrangement of the 
dither matrix throughout the picture, a limit is placed on 
the low frequency content of the output. The technique, but 
for this limit, produces a complete range of frequency and 
tone scale- The most critical artefact that is seen in 
this technique is the arte ficial texture that often accom- 
panies the image. This texture is visible in the dark areas 
of the coat and bow of Lincoln. To develop an algorithm 
which gives both high detail rendition and minimum textural 
visibility is a topic of ongoing research. 

(d) Constrained Average Technique s 

High frequencies (fine details) are obtained by this 
algorithms due to its inherent capability of edge enhancement. 
This detail, though, is at the cost of edge noise. Edges 
appear more blurred than in the case of fixed threshold due to 
the use of a neighbourhood average- 
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(e) Error Diffusion Technique s 

This technique produces the complete tonal scale* The 
spatial frequencies are a function of both the gray level 
being reproduced and the input image detail. There may be 
a spatial translation in the location of edge detail 
depending upon the image content above and to the left of 
the edge. Tho technique produces lower frequencies which 
are more visible than in the case of the ordered dither 
technique. This effect is clearly visible in the pictures 
of Figs. 7b and 7c. 

The control of the spatial frequencies of the intensity 
errors by the ordered dithor technique provides the best 
rendition of the complete tonal scale in the original image. 
This is seen clearly when the dither outputs are compared 
to the outputs of the other techniques. 

Second conclusion reached is that beside the fixed 
level threshold, the constrained Average Technique due to 
its inherent edge enhancing capability will make the detail 
in image consisting of lines and text more visible. 

Jarvis et al [l] have made a comparison of the various 
algorithms and presented it in tho form of a picture showing 
tho comparison of the gray scale capabilities of the various 
techniques. They have reached the conclusion that Ordered 
Dithur bo.idos being a simple technique to implement, also 
gives tho best results, and, perhaps constrained average with 
dither is a strong contestant to being the best. 
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COpiltJIS^THE ffiMSiOW rACTOR.TrilS iS NQM ZKRD 
TMDiSAfES THB CURVINESS .DESIRED, TF ABSCSIGMA)1S '/fiRX 
DARGECF;*G,SO) THE RESULTING ’CURVE IS VERY 
ROySaMAl Up.fHE SIGH Or' SIGMA IHDXCAIES WHETHER 
SiiDRE TflifORHAftOH HAS. BEEN INPUT OR HOT. IF SIGMA IS 
NEGAJXVe THE |nD POINT SliOPES WILE BE DETERMINEO ^ 
iNTERNAyLY.A STANDARD VALUE FOR SIGf<A IS APPROXiMATbr.V 
VALUE. 

ON £ifJTP0f«***»»« ' - 

N/X#Y,5LR1#SLPN AND SIGMA ARE UNALTERED, 

KP AMD IP CONTAIN INFORMATION ABOUT THE CURVATURE OF X 
CURVE AT GIVEN NODES, 

AND ■ ■■ 

S COMTpNS THE POLYGONAL ARC LENGHT OF THE^CURVE, 
subroutine kurvI cSfXrY,sLPi,SLPN,xp,YP, temp, s, sigma J 
INTEGER iN . 

REAL :XCNi,yCN>,||PCNJ,yP(N),TEMP(N),S, sigma 

DEGRAD «I^-*;HS§265/180, 
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RiIht’^SaSd SIDES or trioiagonal linear system 


DY 2 aDELy 2 /OELS 2 ' ' ’ ^ 
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DELS 2 «SaRT((XC 3 )-XC 2 ))«'* 2 -Ky( 3 )-Y( 2 J)T» 2 ) 
DEGS| 2 «fQELSl+OELS 2 

SK»P»¥u ;t?!^*^i 2 Uc 3 »X( 3 ) 

kf Pr 1 'ftifkiy CSY . SXi 

0 lllMl«lQRTtcMN*' 2 >*X(NHll)»* 2 +(y (N- 2 i-/ (NMi ) )** 2 ) 
D|p;iiRfC(iCliMb-XCNh»* 2 icYtNMU-nN))** 2 ) 

OgLp»DgI#MMi+DELM' " 

||jPix(M-l|!c 2 ix?Sil)+Cl*XtN) 


| 7 sy<) 
iibiio 
17700 
17 Ha(’ 
17900 
{ho 00 
IHt 00 
|a2uo 
182*.>C; 
182/5 

IB IDO 

18480 

\85mO 

18 AU{i 

;,87u0 

OBRUO 

189ut) 

190U0 

iuiiio 

192mO 

:,93U0 

.|94U0 

■',9500 

,9600 

',9700 

t9R00 

19900 

20000 

20100 


Performs the mapping of points in the 

A®pSRA"ElllPle(lRESEBflIiloS®4SlB □PPo|e''c|«Po|e3t| 

ARE SPLINE UNDER TENSION AND FUNCTIONS OF THE PQLIGON- , 
AL ARCLENGHT PARAMETER, 

T^'^COMTAiNrrHEAL VALUE Of ABSOLUTE VALUE tjEfS Wg OR 
EQUAL TO l.TO BE MAPPED TO A POINT UN THE 
sK.tv JF T fs iGNOBEU AND THE INTERVAL CO -li, ) IS, NAPPED 
nJiTO THe; ENTIRE CURVE. IF T IS NEGATIVE fftlS INDICATES 
t.iAl the SUBROUTINE H^S BEEN CALLED PREV|OUSL|CNIfH ALL 
other TUPur VARIABLES UNALTERlDj , AND, THAT THIS VALUE OF 
tSu:ESOS the previous, value IN &|SOLU|e VALU|,NI|H 
bcH IMFORMATION the subroutine is AIlE to map THE POINT 
much more RAPXDLT,THUS IF THE OSER Si|P-.l 2 r!J^PM^v 'i« 
SEDUENCE or POINTS on to the same CU|V|/BFFICIENC| is 
GaInId BY ORDERING THE VALUES lNCR|A§|N^^IN, 5 .tAGNI|U 0 E 
and setting the sign of all bUT THB^^FIRST NEGaTIVS 
NCOnIaiNS THE NUMBER Of POINTS WHICH WERE iNfiRPOLATKD 

*°s8PI'*JJrARgLi“cONiAI»llN« tm X-ADD Y^COORDINAfSS 

5uRp|up*IP3«^p^ ^^^'^*** *** "' 




I 






Dir42«oi:t-s*OEr*i 

|}§§5P§ |J?|i!?l2f:;?kps) 


IS)*BCK,J) 


0 f |S‘ t !il?H ' sSaMENf TN IS 

:|H:H 

SU*^»SUM+jSst4S 

IF ABS(T)' IS CREAfER THAH 1., RETURN TERmIWAU 

POINT on-curve 

XS»KtN) , , ^ 

- TSWICU) 

RETURN 

SET UR ANO PERFORM IlfTERPOl.ATlON 


|g«||p(|3fsfJ§Srl|^fl£f]»S|NHD2)/SINHS 
I f C| X £ C I ) lloilil ) «XPt I-l ) ) ♦DEli2 ) - 


i/dels 


UBRQOTtNE PEfFORMS THE MATRIX MULTIPLICATiON 

Si s^MfltliS^TO GIVE c Matrix as the result 


FUPMM'C4Fi6.10) 

Return 

PiviD 


'‘US .>r .irjr ?>!<•. 'iea»!rat.ps cr os s.-i^^tcrrYP j cl-4arf*, 

'I’C i '> r t:t:jOrf* '*") 0 S 5 Hoes to I'll? 

0 j’ '' c 1. 1\' i 1 oy' ;ioyto«r tij'-if'-'* -^oosP 1. ioas ’“uo otiri-ni?! 1:0 * 'yo 

^ 1 > I* ^ 

u ^ :;s' -;» j'j ‘«rtSKf , n j.r-j'iKf 31 j , r? i ) , i 3 • ) 

i ! ■ ’j" a 1 t -i ) , X 1 ( s . , 3 1 ) , / u 0 j , 5 n , i 1 1 ; J , O ’ 

^0 f H J, V, l ;i •■!')■, 0, 0.'! , 1J >,j. j/ 
lobij 'i U'Otl'/t, 7 J 
:!'W,( j '! MO,, f .. ■> 

’’li'it ^ lour? 

';a,oi'"'i'"a. .laca ruoYlao ncir:ill,»i to x-aicis 
0 .1 ) 1=1, 

7 J i j = i , 0 1 

■^ :. 7 f \ M t X U T , J j , K t (, r , ) J , i U t , ’ ) 

':,■) i i'T »w’i. 

' J <0 a 1 , 2 ? 

l>:) to oPtliMi a. 4, HI 

'<('iPOi.Oj 15,01 C 30 -''l.T:.f---'rJOX tM ■ 

-iPJi ^T) 

x( •:p^i^'l')=^’ar i j-oyi ei 

!! JMfT.* Ifi 

Plot a«»?o Hoe as It is rj'"p!t».i 

:a.u. o.'M 7.2, o.?, 

t MfunSMI lo.a OJ 

Z 3 '' iT Jifi, ■ 

SeAOooi' tl,|.or 5 

:j.*osrati, data ruonloa oaraliei to ' 1 -%k^s ■ 

J>1 53 ?ir.,i«^si,,20 

33 43 J'lPjrHTa' , 5 i 

on , Ju ri is-t jo-"?! ( in , 7 , t MT) 

Plfli' 8050 i’ioe as It is co-npotea 

caut<' »b U‘ 3 .jc till, t, j'lrsjp, 2, 0.0, o.2.a.i, 

1 51 ,- 23 . 0 , -.i a “i. 3 , ». 0 ,O, Mftji, V*:RrE 7 V 

C 3 M,l'fiO£ 

Sfi^' i fc'raita 00 fT** fiaa'-e. 

?*Ram£P( i. ”e;‘=’r'^K, '-US^ ) 

5 i 7 P ? '? .1 n 





’tssumeri to 


ttils ilne 


s u,e c t £ j/ t -.'Vt J o a s.l f e ‘ii J • <:• 

r . ■^ Svf. lolt .i.>f.s n£ 

irt:tsr«s at' f‘«^.|allons o 
»■»£, I*; ofics- •.‘/ •‘n ■« 

tHc , t.s.r'- ■•f-ei 

a/.tF, i-j--' =■•!•/ -i-ar^ -"-ro 

cr o »i n t e r - c 1. 0 c n *’4 s b's*,i» fc 

t s • 'i, ao ? 1 1 «‘i I’t* 

»i oj e e 1 1 0 ■ I 0 C . ; t ,0 ? ^ 1 1 ^ 1 ; 
Xo tar^s .0f ro;ra^ttftr||'^ 
rot atari ov' *'■' 
t '.'t I a s a 1 1 a i 'j 


.1 1 

,( I ? ^ 

) r-' . 

^Z 


0 0 / . 

Oi' >i-‘ ■ 
OC 1'.< 

Q.'iS ! 

0 1/ ^ „ 

005 t 
OO'iS. 

oy5f>'‘ 
0’'>5 7 . 

0t'5n-: 
00*5 ‘i.; 

0 0 Cl ' V 

0u6 I 

0Cf.3v 
0 0 6 4 

006 5 
OOodC: 
006 70 

006 8 y 
006^0 

007 00 
0L.710 
007^0 
0073’i) 

>07 40 
>07S:., 

!SI“ 


KKf-' 


V J r 1 3 as 

i ^ i -■( '■);) tat ; , 

D h D r «' V i O M S 


Oil'' 

wO'Jm'SiJ’) 
U'.’Ktjn 
0'.“' I 
On I .' 
0'J‘*4 . 
009J 
009 4’/ 
0096. 
0096 n 
O09 7 m 
O098O 
00990 

0|02i,- 

,v..Ol i) J v.« 

'' p04a 
vOiOSO 
tIOoO 
»J07o 
>|080 
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-V * 0 
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: 1 44 
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,J ‘.J 1 fc/ X t ' « ir u *, 5 I ^ JI jH W t 4; # /% I 9 n 1 « f I, .<• I z \ #» <■ c I* f • . 

1 XSCAbfc, f'-UNti. p^'l, tai-'Tii 

7 ZRij.-% Xt/£\*THf -.ASK, vEHTfc'X ) 

.ias<e'i j'-Oiif-ens tonal plot proaran- vat-'i rotntlnn: 
lals routine ^vtn accepts 3-Dlr.er.s tonal •-'ntn in 
tori-s ns Input f rotate In 3-space fo -any annle 
nnii olot tne prolertiof* ot the resultl'^o luurp 
KV-oldne, Linear interpolotlori Is used uoiuefn 
Tnose lines ot a tlaure wnjcn saoula oe Mriupn 
line are *t\as!<ed« 

rrte T’asitinp tecnniaues usen by ti.is routlcin is onsea on 
t(topte'nlses- 

Ltnes In tne tore-urouna i positive Z-rti r ec 1 1 un > 
are ©lotted oeCore lines In me tensniojon, 

A line or a portion of a Hue Is vasKeri f oMdeo 
It lies ritnin tne reptori noJuded ny oreviousiv 

Sacn call to'f-L0T30 causes one line ot a figure fo oe njoi:t«’'i 
r^o parameters of the plotter ere set on ine Inl. tial coll 

( PI PI ) is tne numrer of plotter increments per loch, 

( i^XPl ) is tfie numoer of liicrements avdllnole accross too 
width pf the paper ( t-dlrection ). . , 

jifnen a new figure is Itiltlated, tne Plotter origin is set 
at tne oottom ot tne pacer nv PLf’riO and snouid not oe 
moved until tne figure Is completed, 

T n 0 H fc ' OS 1? '3 ^ 0 t, 0 S **' 

( 3 is a four dlalt rieclmal liitoqer wnlcn Is used 

to select various input/output options. Tnese diilts in 
decreasing order of mafinlttme, v^lll ne referred to -as v, 

St,. 3f, and z, ■ . , ■ ■ , , ■ , 

I£ V#mS# 0, tne vertlcles of tne current fioure anu tneir 
®tdl*ctlon on to tne i = o plane, will oe stored in a 1b 

,tTf real array l VEHttX ) atid will oe updated as eacn 

„,nt Is Plotted; These co-ordinates are in Incne's ani 
relative to tne current Plotter orJuln, Tne X-X pairs 
‘'■.iRi ’Ordered so that first pair corresponds to the 
'■flfSt point of tfie figure, tne second pair corresponds to, 
th* last point of tne first line, and tne following 
pairs are ordered in a circular fasnlon, roe pairs on the 
y » 3 Plane of the tiyure, then follow in the same order. 

It y = 0, tne VERTEX pararaeter Is Ignored, out snould not he 
deleted, ■ 

ff K = i», the X^components of this line are assumed to be 
exuallv soacpd, and are computed bv 

XC I ) = XOArA(l) -t C I-l )»XSCALK 
fi'neie ( XhA’».A 3 is the initial vaivje in Inches and t XSCALE 3 
is tne soicing oetween points In 'inches. If X.NB.h, the 
x-co'' jotienis of this line are read from an array and 
to fieri nv , ’ 

XI ) a XDrtTA(l3»XSCALE 
4uere XSC'd,iF. Is a scale factor, 

>'•)(' same relations hold for the y-components , 1 

^ ' U 1 3 -a -/DArAd) f (I-D^ySCALE 

and i,t, 1 ,:js,0 

Y( ! 3 a YaAtA(13fySCALE 
ft Z s a, tne Z-components of tnls line are all 
be egual, and ate compated by . , 

Zi I 3. » ZDAfld) f (HLiME-DfZSCAiiK 
.Doere ( 3 Is some ir ' 

’tf Z.siS.O, again we nave 

Z( I 3 « ZDArAtl3*ZSCALE 

When t MLii^E 3 Is equal to one, it indicates the.* Peglnnina 
of a new figure, 4 call to pLOTiDt.with (■ NLINE.H equal to 
4ero .before Initiating a .new figure simulates a line,, drawn 
at the ooctom of. tne oagt* '-..Therefore ofiiy tn.ose portions 
of a line lying .aoove' all .previous lines will^'Cn 
all other parameters are tpnPrea"'Ph such'^a ca'li.' 

( gpi'ixs 3 is the niiffiber of -*.#01111 S'4on tnls line and may- he 
al-tered tro'i' line to line'. .. , 4?. ■ ■ 

f pg’ 3 and 1 i'byXA 1 ar© t'ne tiro- angles (■ in degrees 




P 

a 

qn i j 0 

L 

«P 


tr 

i 

nit 

\ 

al svst 

e 

Qt i 

ixes 

f-clgi f r. .}T ) 


C 0 11 


1 0 

cc 

t S 0 

r3 0 n 

'i¥St 0 


'ts 


1 .A'-'c 1 f e 

(i 

t 0 ^ fi 

: ij 

. 1 * .1 'f t> t » 1 0 <1 


by. 

"i 

arqir 

f 


A 3 

'/ vt -ey, 1 s , 




tats fl 

"■j 'i 

1 system 

•S;;' T-i;o, Piott 


d -f 

\ 

•tj.r,* 




..hi ii-, 5 -.fe o 


i: 

J 

..£ t ..-p 


hO'. 




fr 


he 1 i pu 

fE 

t S 1 

1 rs 

%Z.Tk'.if.PCl0CK 


isa 


'nan ■.!}*. 1 

It 0 


Is, 

t; <?h'. r<^t<hto 


.iV- 


ah"an«.l 

0 

t, ^ ' 

Swwi 

X;.!.--, . C f l'> 


te i 


f 3 n; r e 

Is 

t 'U? 


oh- 


0 t 


e xy-pl 











Itt'r t:.P fpr 

B 



/bee 





f tvi_. 

± 

> 1 * g t’ r gs 

f M fcU y 



K. O 




W1 




( ) and ( ) are t.-te co-orairtatf*s in 1nc:n‘=‘Sf 

relative to the Plotter orlqlt'., to oe useo as tne "ir'-do 
ottneirjdre, ■ ^ 

( XijjlJi’H > is toe ienqth, Ifi Inct^es, to which toe oiot is 
restricted. Any points wMcr. excee^^s this lif’il!', or 
Hfii^s of toe oaoer in the Y-olrection ( -vni ) -ill 
set to to it limit, , , , 

f '^hS'i ) Is an Inteuer array of P + xi/Fo tn ♦f i PI entries ►"ilci- 
Is ase.1 to store toe ■-•a.sk, '(he ceoteots oi l*^1s arr-^y 
sooald oot oe altered ourino the olottlnq of tO/ ^Iven 
P I d 0 r e , ' . 

ili oarateters except ( ,\S^. j and ( VPKitX ) ar™ ret'itno" 

Aet»v98n any t*o calls for tne sa-^e tlouref any naraner.er 
can oe meanindfolly cnanoeri excert ( xiiSnid ), ( '■■’vsk ), an 
( YERiSx J , 

Ii*r£3eR tllGH. UtOHl, ULPLOV; 

XDAUC 1 K iuAlA( 1 ), Zt)AlA( 1 ), ''ASnf 1 ), 

1 t/SRI'SXd) . ., , 

0A?A i.dT, JUtL, SPPT, STHPTA/ -1, -1, -J.tfPiS, -l.ir.f!':'/ 
Ifiiti al i?;atian procedures 

Initialisation procedure tor a nev tliire 
test for special '^Ask u-ooifylnn call 
ir( ^Lioiii.fiu. a ) tiuiu 5oo 
Determine If initial I ?.at ton Is reniUre-i 
ir( ) GOTO 20 

Set plotter parameters 
Pl^£ * lOO^O 
«)1PI * ia90 . 

Heset Plotter orljin to bottom of Pint pecje 
■ I,. «-‘^*i|p!+iuo 

Jftlii. ,!pgotc 0, 0, 3 ) 

Coftpate lengtn of plot paqe In increments 
^ Lim-TX * XUi;Mfd»{»lPl + tO.5 
I'-f C»ItITX t UidITX 

luftlallife masiclnci array over tne lenqtn of tne plot pane 
'.‘viOD 13 K=t r I . ■ 

>v MASX(t<) s i^nr 
COfifin^JS 

l^Ii’ * -1 ■ ■ 

Set ttie necessary Indicators for tne first line of a nejK 
f ioare 

l.'d s -1? I a O'' - ». 

rnput tyoe and v'ERfSK Initialization 
)eternlne If Initialization is required 

fF( j*/<yz,£g. ivxxz ) goto 7o 

Set Indicators for types of input data and savlno vertlcles 
NKU = IVKVZ 

rOZ = IJ IDDI a it IMDX a If II^DV = I 
TFt JVXfZ, uT, 1000 ) GOTO 30 
t 'i di/ = 2 

r/nz = jvxxz - looo 

rf( Ji<XyZ,!jT,iy0 ) GOTO 40 
T.‘ OK - 2 

Tl/XVZ = lYXXZ - lOD 

U'C lyXfZ.GT.lO ) GOTO 50 
i;i3X S 2 

JtfXXZ * JVXl'Z - 10 ■ ■ ■ ^ , 

IF( J¥XtZ,GT,l ) GOTO 60 

INOZ * 2 " ■ ■ . : v ■•'iFF.r. , A'KtF 

jni’Z » Tvxtz .. .. ' 

Rotation lnitlaiii$3tlQn • 

Determine If inltiallaatlpn ,is require a 
IFC PHI, EQtSPHl.ASlD. THETA. CO.STHRTA 3 GOTO »D 

=°piii%liiPs“8 

CTHETA^« COSFC :0,0!74&3a$25#TMEfA' ■ 

Ml .= CP.JI . '.:;Fv,.r 


)i 3 > ' 
u i 3 I’; 
Jl U I 
n I 3 n 
;n ’ “f •'■ 
dll' „• 
■ji If 
Ui U-i 
■j ! 4 3 
U 1 d ■( / 
u I I n ^ * 
■j i 4 h a 
v) 1 '+ 7 
w I -I H ■' 
D I ID 
ilSht. 
D l P 1 

ul 5i^ 
Old 'id 
O-I 55D 
015 6 v 
yi5?u 
OISHO 
01590 
0J6D0 
01 610 
01620 

Ollll 

III 


.»17RD 
J I 7 9 . 
U I ' * 

0 J ri 1 . 

0 1 n / j 
013 3 , 
OIS'K* 
aiHbu 
01860 


1910 

i92u 

1930' 







Sf a XDATAtD + XSCAbE, 
jjj’j i nOf 13D ) , n-iy 
^ a YDArA(l) f (T-1 
”jr3 I 43 

‘i a vaATAUi + iSCAl.K 

33i‘.J (. 150f 150 ), 1^07 

?. ■= Zi)ATA(iJ t (.^LIwh:-! j^ZSCAL?' 

3JT3 1?U 

Z a ZOATAtD^ZSCALe 
idta rotation 

Hit a AllAX + AU»Z + X.BPF 
X< a XXX 

IX a t-nxC XXFPIPl + 0,5 ) 

{It = A21FX ♦ A23FZ f YREF 

Vf a {it * A22FX 

II a IFIXC Xr^PiPT + 0.5 ) 

Pf»5trlct Cl^are to plot oaae 
U'l rx.bS.J ) IX si 
IF( IX.GI.LI-il'fX ) IX a L|: lTx 
ifC Tx.ur.ii ) IV = i) 

£F( tV.GT.J-llPT 1 lY * OYPi 
IF( k.ne.i ) Goro 2r.o 

( b'Ow ) Is tne position of tne cre^/lnus point «1tn rt'sooc 
me Aasic 

tl aoove the ftasK 

0 'Titnln the llfalts of tne nasF 

-I oiiow tne :r«as<<. 

Procedure tor Initial ooint of each Line 

bocate initial point sslth respect to tne inasK tnen 
updatt me masK 
. : tbDsi a ix t IX 
. * bOF - 1 

„ , *|ba«| a ‘lASKC LO*« I 

n HHiSH a 'lASKl HIGH ) 

■,;r|F( MrllGH-IY ) 200, 210, L 8 .) 

■'■IFC MbOH’-IY ) 100, 230, 220 
"4o:Gbo * 0 
. SJT3 240 

" 11 ASX C HIGH ) a ■ lY " 

^ . IFC , MbO»v.eQ.-i ) maskl low.) a lY 

b 3C ObO a I , . . 

G3rn 240 

iAS*f( mOO ) a IY 
bOCObO = -1 

iova tne raised pen to the initial point 
: Vub r, >!,:;! i u, u, 3 ) 

ir a u; I¥ = lYl HHEF = iY 

''core vsrtlcies it requested 
(FL C'rll.Kj.l ) GjTU 530 

I 't 0 E X a f „ G i f 5 

'/iP IFXC J '.DCX ) a • 

i'Exc i yoEx+i ) = vv ■: 

I-:bl'fiX(l'vU£X + d| a XXX 
lERTSXClMOex+O) a VYY 

IFC Hbr^E.'^F.t J GOTO 530 ' 

yERl'SXC I ) a XX ■' ■ , : 

OERTSKC 2 a VV 

VEBTSXC 9 ) a XXX "OR m 

YSRiexC 10 J a^YVY ■ , 

GO IQ 530 ' . , 

Saectal case wnere cnanqe In X co-ordinate , IS 
A' special orovlslon Is tr.ade at this point sp t|\at a line 
Aas<' itself as long as tae K coordinate resnains 
constant , . . , , , ' 

IFC IX.nE.JX ) GQX.O 260 ' ■■ ^..R' 

JV •■* if . '■ 'R .,RR'-- .■ • ,v.: ." , 
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0X3 
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0 2 3 3 . 
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02'U:„ 
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02410 
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0329.; 
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033 3c 
0334) 
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0 33OW 
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03390 
03400 
03410 
03420, 
03430 
03440 


tiiir'Ji. ^uxiit i& ■«[ 1.1,11*1' i UTT i'ia» 


i r % u w-' ^ j u j "♦ 'j f a I ' f a .a 

■fa© rurrent. ooiat Is belov; the mds*c: 

f ^ Q s " ^ I 

i>'( liJCOf.O ) 3Ja, 4Su, Ur., 

Plotfraii aboye tne masK to within toe ’tasi*; 

.lf( '1HH’3»3.LL,1YREK ) Qkhh TfLO'^f JX, .liUGti, 

'3jr3 350 

olol oelow tne "^asK to witnfr? tae 'nasK. 

IfC 4L0.3.r,£,I!t«fc:F ) CaLI, IPiiOTt Ja, 'iL.nw, 2 
Plot fton within the 'nask to witrin tne mask 
car.b IPbon JX, Ji, 3 ) 

Giro 5'2D 

Plot trom DBlow the -nastc to aoove ti\e 'nask 
tF( 4uO^-IYREF ) 37a, 3t.a, IRf 
P.lot fro'n wltnlh the mask tu above ih.^ u.iSk 
rf{ 4 HI an- 1 '/REF ) ivi, , lya 

GAGb rPbOTC JX, 41, ;n-^, ? i 
GALb tPLOK Ja, 'iHiGa, 3 ) 

G3I‘3 430 

TF( MrtIGH,gQ,-l ) GOri. 43U 
3b3m .* HIGH - 2»IHCX 
TF( HASA(ai.l)Hl)-,Ty ) '*2!/, 420, 410 
CAUb IPliOTI JX, JY, 3 ) 

GDTD 430 

IPGOTC JX-I<«CX, HASK(ni,nHX:) , 3 ) 

Plot from aoove the ;fiasK; to aoove tne masK 
■’ HkSKi HIGH 3 S JY 

IF( ) XASK{ LOW ) * jy 

iaroslo''" JX" ^ ’ 

PUW FROM ABOYE fHE HASK TU BE!. 04 I HE HASK. 

» TF^ HHISH-IYREF ) 450, <ioO, 460 

■ f rom , within ■ the mask to-uelow the nask 
HtOil*iyREF ) 470, 47u, 4R0 
.■■:yZkhh XPLOTC JK, v-^HIGH, 2 ) ■ 

i’^3L0w a - ^’f'PICX ^ 

irc 4A3iCt;jr,0{,a43-jy ) 490, 500, 50D 
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